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Abstract. We discuss the behaviour of a thin soap film facing a frame element: the 
pressure in the Plateau border around the frame is lower than the film pressure, and the film 
thins out over a certain distance X(t), due to the formation of a well-localized pinched region 
of thickness h(t) and extension w(t). We construct a hydrodynamic theory for this thinning 
process, assuming a constant surface tension: Marangoni effects are probably important only 
at late stages, where instabilitites set in. We find X(t) ~ t 1 ^, and for the pinch dimensions, 
h(t) ~ t~ 1/2 and w(t) ~ t" 1/4 . These results may play a useful role for the discussion of later 
instabilitites leading to a global film thinning and drainage, as first discussed by K. Mysels 
under the name "marginal regeneration" . 



Early experiments by K. Mysels and coworkers [jjj showed that a vertical soap film, sus- 
pended on a frame, (and made with "mobile" surfactant) thins out by nucleation and growth of 
black, thin spots near the Plateau borders. They called this process "marginal regeneration". 

There are in fact (at least) two steps in marginal regeneration: a) The pressure in the 
Plateau border is lower than the pressure in the film. This thins out the film near the border, 
and leads to a "pinch" . b) The pinched state must have an intrinsic instability leading to 
the black spots. The two steps are very different: the pinch can occur at constant surface 
tension — i.e. without any Marangoni effect. On the other hand, the later instabilities are 
probably triggered by Marangoni flows, as pointed out by a number of authors 

Our aim in the present note is restricted to the first step, i.e. the description of the pinched 
state with its dynamics. Pinching has already been studied in connection with the elimination 
of dimples in the coalescence of drops HQ or in the drainage of thin films |§,[)). However, in 
these problems, the extension of the dimpled zone is prescribed and fixes one of the spatial 
scales involved. Our case is different: we consider a semi-infinite film (of initial thickness eo) 
facing a straight Plateau border. At t = Q the film begins to pinch, and is perturbed over a 
certain distance A(t) increasing with time (as we shall see A(t) ~ i 1 / 4 ). Also the width w(T) 
of the small pinched region decreases with time (w(t) ~ i -1 / 4 ) and so does the thickness of 
the pinch h(t) ~ t~ x / 2 . It seems important to know these scales in detail before embarking 
into the second step: all data on dimples suggest that instabilities can only occur after a long 
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Fig. 1 - (a) Soap film held in a vertical frame, and xy cross-section plane, (b) Sketch of the initial 
shape in the cross-section plane (here, a case where the liquid totally wets the frame material). As 
the film in the central region is very thin (uniform half-thickness eo < r), it has been represented 
simply as a line. 



stage of pinching. This is why we study the pinch in the simplest condition of a constant 
surface tension. In this CclSC, clS will be shown in this article, a complete, explicit solution 
involving the the two scales A and w can be constructed. 

We thus consider a soap film, held in a vertical rectangular frame as in fig. |l|-a, and we wish 
to describe the time evolution of the film profile in a horizontal cross-section, neglecting in 
the process all drainage in the vertical direction (vertical drainage in "mobile" films is tightly 
associated with the second step of marginal regeneration, and will intervene significantly only 
once the instabilities have risen up). The cross-section profile is therefore treated as a one- 
dimensional thin film of half-thickness e(x, t) at position x and time t. 

The initial shape of the film has the following general features |jj (see fig. |l]-b) : joining the 
frame legs is the "Plateau border" of size r, where the thickness is of macroscopic dimensions. 
In this region, the curvature of the film is prescribed to a fixed value c. Further away from 
the frame legs, there is a transition region, where the thickness decreases steeply. And, finally, 
we reach the zero-curvature central region, of constant thickness 2eo, which extends over 
the major part of the frame width. The system displays in addition a natural separation of 
spatial scales. Firstly, the width of the frame (e.g., 10cm) is supposed to be much larger 
than the Plateau border extension r (e.g., 1 mm). We shall thus assume that the width of the 
frame can safely be taken as infnite in regard of the length r, leaving us with a semi-infinite 
film. Secondly, the Plateau extension r itself is much larger than the limit thickness eo in 
the central region (e. g., 10 /im). Our study of the film dynamics will therefore rely on an 



asymptotic approach valid in the limit e <S! 1, where e represents the aspect ratio of the 
film: 

e = e /r. (1) 

The evolution equations are classically written in the lubrication approximation. The axis 
are shown on fig. [lj-b, and, for simplicity, we choose the origin of the a;-axis at the position 



of minimum thickness, in the pinch ( ). Assuming that velocity of the liquid vanishes at the 



( 1 ) This kind of approach has been initiated and developed by Jones and Wilson ]7| in the resolution of a 
closely related coalescence problem. 

( 2 ) We assume that, once the pinch is formed, its spatial position does not change anymore, or, at least, 
that its motion remains of subdominant importance to the solution. This hypothesis is indeed supported by 
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film/surfactant layer interface, the flow is of Poiseuille type and yields a parabolic velocity 
profile. With a Laplace- Young pressure p = p a ; r — je xx in the fluid, we find that the flux Q 
in the ir-direction is given by Q = ^V*e 3 e xxx . Here, V* = 7/77 is a characteristic capillary 
velocity and is defined as the ratio of the surface tension 7 by the viscosity of the liquid rj. 
(7 might depend on the vertical elevation of the considered cross-section, but is otherwise 
constant in the xy-plane.) Inserting into the continuity relation 2et + Q x = , we find the 
governing equation of the film profile: 

V* 

e* + — (e e xxx ) x =0. (2) 

Note that microscopic interactions Jl(J like van der Waals, or electrical interactions, have 
not been included in eq. (g), since we are mainly interested in thicknesses larger than their 
range of action. We will simply consider that once the film thickness reaches small enough 
values (~ 100 A), these additional interactions are globally repulsive and prevent any further 
thinning. 

It is more convenient to use the reduced variables (denoted with capital letters) X — x/r , 
E = e/eo , and T = t/t le \ ax , where t re iax = y^ rA l e o is a typical relaxation time needed for 
a bump of height eo and width r to be leveled off by capillary flow P)[ In these reduced 
variables, eq. (^|) takes the dimensionless form 

E T + (E 3 Exxx)x = 0. (3) 

From the initial shape described above, a pinch has first to form. Describing this process 
requires the precise knowledge of the initial state. However, it is expected that the formation 
stage is quite rapid on the time scale of i re iax, and, moreover, the details of it appear in fact to 
be unessential to the longer term dynamics of the film. We may nevertheless understand why 
only a localized region thins out at the beginning and starts to pinch: within the initial profile, 
the pressure gradient is non-vanishing only over a very restricted zone, inside the transition 
region (there, the curvature decreases from the border value c to the central value zero). 

We will hereafter simply assume that this early stage leads to the formation of a pinch 
precursor of reduced thickness. Without the need of further knowledge, we are then in a 
position to describe the evolution of the soap film. In order to do so, we solve the equation 
governing the dynamics (eq. (||)) as follows: we divide the film into three spatial regions 
(inherited from the three regions in the initial profile), and look for solutions in each of 
them, imposing in addition that these solutions all match together (in the sense of symptotic 
matching). Starting from the frame leg, the three regions are: the Plateau border, (at negative 
values of X), then the pinch region that lies in the vicinity of X = 0, and, eventually, the 
central region, which extends over the positive values of X . 

The profile in the Plateau border is easily determined to leading order in the asymptotic 
parameter e: as the border contains a great mass of liquid compared to the rest of the film, it 
remains unchanged during the whole evolution. In particular, its curvature remains the same 
as initially (i.e., e xx = c), maintaining thereby a constant pressure drop as compared to the 
pressure in the central region. This feature will naturally prove essential to the film dynamics. 

For clarity's sake, it turns out being more convenient to pursue with the description of the 
central region (X > 0). The existence of a pinch at X = 0, with a sudden thickness decrease, 

previous work on dimple formation in analogous situations ^||s|]. 

( 3 ) Here is a simple way to see this: the capillary flux, as used in establishing eq. (j^J), is Q ~ V*e 3 e X xx ~ 
V*e 4 /r 3 . Then, i ro lax appears as the time needed to empty the volume of the bump with such a flux, i.e., 
*rdax — e r/Q. 
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involves that the central region, initially flat, is pushed out of equilibrium, and that a flow 
must take place in order to relax, to the greatest possible extent, the curvature induced by 
such pinching. But this flow cannot easily pour liquid through the pinch region, whose small 
dimensions strongly reinforce viscous resistance. Thus, the relaxation flow is mostly directed 
towards the far-field zone "at infinity" (X ^> 1). In a more formal way, the relaxation flow 
Qrciax towards the far-field zone is of order unity, whereas the flow poured through the pinch 
Qpinch is only of order e (as will be shown later on). Thus, in the spirit of boundary-layer 
techniques, we look inside the central region for a solution at zeroth-order in e (a so-called 
"outer" solution), that is to say that the relaxation process is solved with a boundary condition 
of zero flux at J < 1 (in the vicinity of the pinch region). The complete set of boundary 
conditions for the central region solution writes: 

(i) E\ x<<l =Q 1 (ii) E 3 E XXX \ X<<1 =0, (iii) E\ X>>1 ^1. (4) 

The first condition imposes that the film thickness must vanish (at zeroth-order in e) at the 
approach of the pinch, the second one is the zero flux condition, and the third one ensures 
that the far-field film remains unperturbed and retains the initial thickness eg. 

We propose a solution to eqs. (||) and (||) in a self-similar form, with the self-similar variable 
U = -j^X/T 1 ^. Setting E(X,T) = F(U), and inserting into eq. (§), we obtain the following 
equation for the self-similar function F: 

UF' = {F 3 F'"y . (5) 

A study of the film profile in the pinch region (see further down) shows that a matching with 
the central region is possible if the behaviour of F near the origin is either linear, parabolic, 
or power-law (with exponent 3/4). However, it turns out that only the linear behaviour is 
compatible with eqs. (|J) and (JsJ) , i.e., F ~ AU for U <C 1, with A a constant. The complete 
function F can be constructed numerically with the help of a two-sided shooting procedure 
based on analytical expansions at U -C 1 and U ^> f (we omit technical details for conciseness). 
The resulting plot is given on fig. ^|-a: as one can see, F increases from zero to unity over 
a range of a few units in the variable U, and then presents decreasing capillary oscillations. 
The shooting procedure also determines the numerical value of the slope of F at the origin: 
A~ 1.591. 

The important physical consequence of these results is that, because of the presence of the 
pinch at X — 0, the film in the central region is also thinned out, and retrieves the far-field 
thickness eo over a distance X(t) ~ rT 1 / 4 ~ i 1 / 4 , which increases with time. 

We are thus in possession of a self-similar solution of the problem in the central region. 
One may though wonder about the relevance of self-similarity in describing the real physical 
solution. Actually, if we take a version of eq. @ that is linearized around E = 1, and solve it 
analytically with the Green function method with boundary conditions similar to eq. (|J), we 
rigorously find that the solution indeed takes the above self-similar form. We conjecture here 
that this feature can be safely carried over to the non-linear case. 

We next consider the solution in the pinch region. As stated earlier, once the pinch is 
well-formed, the suction from the Plateau border is able to drain liquid through the pinch 
out of the central region, but only with an asymptotically small current Q p i nc h- However, the 
dimensions of the pinch are so small, that we can reasonably take this current to be spatially 
uniform over it ||0,|9| (but not constant in time): Qpinch = Qpmch(T')- Therefore, to leading 
order, the profile of the film in the pinch region obeys the following uniform current equation, 
rather than the full eq. (||): 

E 3 E XXX = Q pinch (T) . (6) 
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Here again, we try a self-similar solution of the form E(X, T) = h(T)S(rf) , with 77 — X/w(T) . 
Such a solution has been proposed in references [Q,^). In this solution, the functions h(T) and 
w(T) represent, respectively, the typical height and width of the pinch, and will be determined 
by matching with the neighbouring Plateau and central regions. Replacing the self-similar 
form in eq. (§), we get S 3 S"' h{T) 4 /w{T) 3 = Q pinch (T) , from what we deduce that S 3 S"' 
cannot bear any dependence on 77 (this would bring in the l.h.s of the equation a dependence on 
X that could not be balanced by the r.h.s.). Thus, the differential equation on S is necessarily 

S 3 S"' = -a , (7) 

where a is a positive constant (the sign has been chosen in view that Qpinch , oriented to- 
wards the negative coordinates, is negative). At this point, a is still unknown, but shall be 
determined later on. 

Let us now apply the matching requirements that the self-similar solution must fulfill, 
starting by the match with the Plateau border. A study of the asymptotic behaviours permit- 
ted to S by eq. (f7|) for large 77 (either positive or negative), opens three possibilities, namely 
linear, parabolic or power-law with exponent 3/4. The Plateau border imposes a constant 
curvature e xx = c, or in terms of the reduced variables, ExxIx^q^ = rc / £ - Thus, on the 
pinch side, S must be asymptotically parabolic. Equating the Plateau curvature with the 
pinch curvature Exxl^^loc = ^ m S" h(T)/w(T) 2 yields moreover the two conditions: 

(i) h(T)/w(T) 2 = rc/e, (ii) lim S"'(r?) = 1. (8) 

77 — *— 00 

We proceed further and match now with the central region. As stated above, there are three 
possibilities of asymptotic behaviour of S at the approach of the central region, but on the 
central region side, we found that only the linear one is consistent for the function F with 
the boundary conditions and the differential equation (|j), giving E(X, T)|^^ r 1 al ~ AU = 
■^X/T 1 / 4 . In the pinch, S(rj) must accordingly have an asymptotically linear behaviour (for 

7/ ^> 1), which leads to the corresponding expression of E: E(X,T)\ 1 ^ 1 ^ 1 = h(T) r\ lim S' = 
X h(T)/w(T) limS". Matching these central and pinch behaviours together, we find that two 
new conditions must be met: 

(i) h(T)/w(T) = ^=T- 1 /\ (ii) lim S> (77) = 1. (9) 

With eqs. (§-i) and (§-i), we are now able to extract the expressions of h and w: 

KT)=e^-T-y\ t B (T)= e -A-T-V4. (10) 

We immediately notice that, as announced before, h and w, the typical dimensions of the 
pinch, are of order e and hence asymptotically smaller than the dimensions (of order unity) 
that prevail in the central region. Furthermore, we see that, as time passes, the pinch becomes 
thinner and smaller in extent. We may also note that the time dependences in eq. (|l0|) are, 
not too surprisingly, equal to the ones found in closely related drainage problems || [7], [J . 

To complete our solution of the profile in the pinch, we still have to exhibit a function S 
that has the correct shape, as dictated by eqs. (0), (||-ii) and (^-ii). Such a function has indeed 
been computed previously At large 77, S admits an expansion around the imposed linear 

profile of slope unity, from which we have started numerical integrations towards negative 
values of 77. The integrations have shown that complying with the requirement that S" equals 
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Fig. 2 - (a) The self-similar function F(U) describing the film profile in the central region, (b) The 
self-similar function S(rj) describing the pinch. 



unity for large negative rj ascribes to the parameter a in eq. ([7|) the value a ~ 1.21 . The 
minimum value of S (at ij = 0) is then found to be S m i n ~ 1.52 . A plot of S is provided in 
%• |b. 

With the knowledge of the self-similar functions F and S, we have now at disposal a 
complete description of the profile of the soap film, once the pinch is formed (and under 
the basic hypothesis that it remains later on fixed in position). We may at this point give 
justifications to some assumptions that were made in the course of the resolution. 

First, it was assumed that to zeroth order in e, one was entitled to solve the central region 
profile with zero flux towards the pinch (eq. This can be justified as follows. In the 

central region, the current generated by the relaxation process is, by the very definition of 
the reduced variables, of order one: Qreiax = 0(1). On the other hand, the expression of the 
current in the pinch was given when establishing eq. (Q): Q P i nc h = S 3 S"' h(T) A / 'w(T) 3 . Since 
both h and w are of order e, we deduce that Qpinch = 0(e) <§C Qrciax , and is thus negligible 
to leading order in the central region. 

A second approximation was that, since the pinch was of small dimensions, the current 
Qpinch could be taken uniform inside it. Indeed, using the continuity equation, one can 
estimate the thereby neglected spatial variation <5<3 p i nc h as the variation in time of the volume 
enclosed in the pinch region: we have <5Q p i nc h — gj^M = C(s 2 ) ^ Qpinch , as had been 
assumed. 

In addition, it has been assumed, when establishing eq. ^ , that the velocity field goes to 
zero at the film/surfactant layer interface, but also that the surface tension 7 keeps a constant 
value throughout the film. This cannot be strictly true: For the surfactant layer to be at 
rest, it must develop a surface tension gradient to counter the friction force generated by the 
underneath liquid flow. We can estimate the required gradient by equating it with the viscous 
stress at the surface: (1/70)07 / dx — {T]/lo)dv/ dy\ y=e = ee xxx (the equations were divided 
through by 70, the constant tension that would be obtained in the absence of liquid flow). 
Using this formula with our solution for e(x, t) , we find that the relative variation of the surface 
tension across the pinch is ^7/7o|p inch = JpinchCVTo^T /dxdx ~ u>(e 2 /r 2 )(h 2 /w 3 ) = 0(e 2 ) . 
Similarly, we find that the variation over the characteristic length r in the central region writes 
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&lllo I central = ^( e2 ) • Thus, the gradients required to cancel the flow velovity are indeed very 
small, and the surface tension can safely be taken as constant (at least, in the study of the 
leading-order dynamics of the film). However, we should point out that this result is by no 
means contradictory with the later destabilization by Marangoni flows, as these occur in the 
vertical direction. 

A final remark is in order concerning the late stage of the film pinching. As mentioned 
earlier, we expect that, at small scales, repulsive microscopic forces stop the thinning (with, de- 
pending on the electrolyte concentration in the film, different possible scenarii || ) . The corre- 
sponding time t s top can be roughly estimated, if we say that the stop occurs for e around 100 A, 
and use eq. ( |T(]| ) with typical numerical values (for instance, eo — 10 /im, r ~ 1/c ~ 1mm, 
V* ~ 10 2 m/s): we find t stop ~ 10 3 s. We should caution that, in regard with experiments 
on marginal regeneration, this is a significative duration, for which one would then probably 
have to include the study of the second step of marginal regeneration (Marangoni instabili- 
ties). The thinning process might therefore not be allowed to proceed by itself as far as i s top- 
We can also estimate the distance over which the central region has a thickness departing from 
eo (assuming that the thinning keeps on until i s top): A s top — Ktstop/irelax) 1 / 4 — 3r ~ 3 mm. 
This distance is much smaller than the frame width, justifying our description of the film as 
semi-infinite towards the center of the frame. However, in situations where the frame width is 
not so large, our central region solution should be regarded as transient: when X(t) becomes 
of the order of the frame width, we expect to retrieve the solutions in bounded geometries 
that have been described by several authors (for instance, a quasi-static parabola in 

the central region). 

* * * 

We are pleased to thank H. A. Stone for very valuable help and A. L. Bcrtozzi for instructive 
discussions. 
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